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ABSTRACT 


The r.m.s. field stren^hs of electrostatic and electromagnetic 
turbulence in the earth's bow shock, measured in the frequency range 
20 Hz to 200 kHz mth IMP-6 satellite, are found to correlate with 
specific solar vdnd parameters measured upstream of the bow shock. 

The largest r.m.s. field strengths of electrostatic turbulence (200 Hz - 
4 kHz) occur when the upstream electron to proton temperature ratio 

is large, and when the proton temperature T^ is small, indicating 
that the mechanism for generating electrostatic turbulence in the bow 
shock is more efficient when lower upstream proton temperatures occur. 

Wo substantial correlation is found, between the r.m.s. field, strengths 
of electrostatic turbulence and three upstream parameters commonly 
used to classify the magnetohydrodynamic structure of the turbulent 
bow shock: the Alfven Mach number the ratio of particle pressure to 

magnetic field pressure B, and the shock normal angle i|r(B,n). The 
strong correlation with T^/T^ and T^, and the lack of strong correla- 
tion with M^, p, and ijf(5,n) indicates that the strength of electro- 
static turbulence in the bow shock is determined by the kinetic 
properties of the solar wind plasma rather than by its fluid proper- 
ties- ; 

The largest r.m.s. field strengths for electromagnetic tur- 
bu.lence (20 Hz - 4 kHz) occur when the upstream particle density W 
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IS large and when the shock normal angle i^(B,n) is closer to 90% 
supporting a previous conclusion that whistler waves comprise the 
electromagnetic tm-bulence in the bow shock. Electric fiald turbu- 
lence, composed of both electrostatic and electromagnetic fluctua- 
tions, correlates with the upstream parameters T^/t , T , and tjr(^,n) 

in such a way as to imp]^ that mode coupling occurs between electro- 
static and electromagnetic waves . 

A bread spectrum of high frequency (3 - 30 kHz) electrostatic 
turbulence typically observed in the leading edge of the bow shock 
IS interpreted as indicating the region of electron heating. Deeper 
within the shock transition the intensity of low frequency (< 3 kHz) 
electrostatic turbulence greatly increases to form a broad peak, 
centered between 200 - 800 Hz, and is interpreted as corresponding 
to the region of maximum proton heating. The characteristic develop- 
ment of the electric field spectrum through the shock transition 

indicates that strong coupling exists between the electron and proton 
heating processes. 


k 


I . irJTRODUCTION 


Several studies of the earth's bow shock have shown that a 
complex structvtre of wave-particle interactions govern the heating of 
solar wind electrons and ions as they stream across the shock transi- 
tion. Fredricks et aj.. [I968, 1970a, b] have shown that electrostatic 
turbulence at I.3 IcHz and the scattering of protons correlates iri.th 
the gradients in the magnetic field at the shock front, thus 
indicating the presence of a current -driven instability- Montgomery 
et ^ [1970] and Forraisano and Hedgecock [1973a, b] have shown that 
electron thermalization occurs in a thin region upstream of the bow 
shock, followed by ion thermalization in the main shock transition. 

Solar wind ions are reported by Neugebauer [1970] to be substantially 
decelerated near the upstream side of the bow shock, possibly by a 
charge separation electric field. Neugebauer et al. [I97I] have 
reported the observation of EIF magnetic field oscillations (10 - 1000 Hz) 
correlated with superthermal electrons ( > 100 ev) in the shock mag- 
netic field gradient. Holzer eb ^ [I966, 1972] and Olson et ^ 

[1969] have discussed the spectrum of magnetic field fluctuations in 
the bow shock below about 100 Hz which shows a broad noise spectrum 
with a peak near 3 Hz, and includes waves that propagate in the shock 
frame both upstream and downstream. Fairfield [197k] has studied 
whistler waves in the frequency range 0.5 - U.O Hz which are detected 
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ahead of the how shock and may result from resonant interactions with 
upstream electrons. The spectrum of plasma wave turbulence in the how 
shock, as measured mth the University of Iowa plasma wave spectrum 
analyzers on the IMP-6 satellite, was reported in Eodriguez and Gurnett 
[1975]* ^-Jhistler waves and electrostatic waves are shown by Rodriguez 
and Gurnett to be the major components of the bow shock spectrum be- 
tween 20 Hz and 200 kHz. 

In our previous report [Rodriguez and Gurnett, I975], herea^er 
referred to as paper I, it was shown that the characteristic electric 
field spectrum of the bow shock can be resolved into a low frequency 
electromagnetic component (20 Hz - 200 Hz) which decreases monotonicaliy 
with increasing frequency approximately as ±0*5) higher 

frequency electrostatic component (200 Hz - 4 kHz) associated with a 
peak in the electric field spectrum between about 200 Hz to 800 Hz. 

The characteristic magnetic field spectrum of the bow shock was shown 
to decrease monotonically with frequency approximately as ±0*5) 

and to display an upper cutoff frequency near the electron gyrofrequency. 
By taking the ratio zj Sg of simultaneously measured electric and mag- 
netic field energy densities in the bow shock, it was shown in paper I 
that Sg/sg increased nearly monotonically from values near 10“^ - lo"^ 
at 20 Hz to 10 - 10 at 1 kHz. The values of Sg/sg (proportional 

to n where n is the index of refraction) at low frequencies and the 
observed upper cutoff frequency near the electron gyrofrequency in 

the shock magnetic field spectrum led to the conclusion in paper I that 
in the range 20 Hz - 200 Hz the electric and magnetic field 

turbulence in the shock is caused by electroiragnotic whistler waves. 
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It was also concluded that the large values of sj/sg at frequencies 
greater than about 200 Hz showed that these hi^er frequency waves are 
almost completely electrostatic and. are associated with the peak 
centered between 200 - 800 Hz in the shock electric field spectrum. 


The frequency range (20 Hz ~ 200 kHz) of the shock plasma mve 
spectra discussed in paper I includes all the characteristic plasma 
frequencies for electrons and protons except for the proton gyrofre- 
quency (~ 0.1 Hz). The heating of solar vdnd electrons and protons 
in the bow shock must result in the self-consistent generation of a 
spectrum of turbulent electric fields. It is expected that such a 
spectrum will be broad enough to include most of the characteristic 
plasma frequencies since these are the elementary excitations through 
which the particle velocity distributions can be broadened. There- 
fore, it is assumed that the electric field spectra of paper I and 
of this study are closely related to the dissipative processes that 
occur in the bow shock. 
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II. IMTEKSITY VARIABIIiITy OF SHOCK TURBULENCE 


A. Dynamic Range 

The variability of plasma wave turbulence in the bow shock is 
indicated in Figures 1 and 2 (similar to Figures 7 and 8 of paper I). 

In Figure 1, two kinds ol’ electric field spectra are shovm, spectra 
measured with a 5«12 seconds averaging time (averages) and spectra 
measured with a 0.1 second averaging time (peaks). One measurement 
of the peak spectrum is obtained within the time period of one measure- 
ment of the average spectrum. (A description of the University of Iowa 
plasma wave experiment and the spectrum analyzer used to measure the 
electric and magnetic field spectra in the bow shock is given in 
paper I.) The left hand side of Figure 1 is an overlay of average 
electric field spectra that were measured in 36 crossings of the bow 
shock; the right hand side is an overlay of the corresponding peak 
electric field spectra. Figure 1 thus illustrates the dynamic range 
of intensities that have been sampled for electric field turbulence in 
the bow shock. For example, on different shock crossings the average 

electric field spectral densities at 1 kHz ranges from a minimum of 

2 -2 -X -"Q 2 -2 “X 

3 X 10 volts m Hz" to a maximum of 6 x 10~^ volts m~ Hz , 

a dynamic range of over four orders of magnitude. A similar range of 

variation is seen in the peak spectral densities. For each spectrum 

of Figure 1, a total r.m.s. field strength E 

° rms 


may be computed by 
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integrating the spectrum from 20 Hz to 200 kHz. The range of E so 

rms 

obtained is (4.0 - O.19) x 10“‘ volts for the average spectra and 

(23. - 2.1.^ X 10 3 volts m ^ for the peak spectra. The electrostatic 
and electromagnetic components of the shock spectrum which were des- 
cribed above are evident in the characteristic shape of the spectra 
of Figure 1, with the best resolution of the two components occurring 
in the average spectra of intermediate intensity. For average spectra 
of greatest intensity it can be seen that the peak centered between 
200 Hz and 800 Hz broadens out enough to nearly fill in the minimum 
in the curve at about 200 Hz where the two components of the spectrum 
meet. For the average spectra of smallest intensity, the peak nearly 
disappears . 

Figure 2 shows the average and peak magnetic field spectra 
corresponding to, and in the same format as, the average and peak 
electric field spectra of Figure 1. There is evidently a smaller 
dynamic range for the magnetic field spectral densities (about two orders 
of magnitude at 100 Hz), and except for a few peak spectra, the shape 
of the spectrum has the characteristic monotonic decrease with frequency 
associated with th i component below about 200 Hz in the electric field 
spectrum. The plasma wave turbulence associated with the broad peak 
in the average spectra of Figure 1 is clearly electrostatic turbulence 

since there is no corresponding peak in the magnetic field spectra of 
Figure 2. 


Correlation Ihrameters 

Since plasma shock waves are usually categorized in terms of 
such upstream parameters as the Alfven Mach number M^, the ratio of 
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particle pressure to magnetic field pressure 0, the ratio of electron 
to proton temperatures T^/T^, and the shock normal angle iJ((B,n), it 
is of interest to investigate the relationship between these parameters 
and the intensity variability of bow shock turbulence indicated by 
Figures 1 and 2. The upstream parameterr, used in this study are 
derived from the measurements of two other IMP-6 experiments, the Los 
Alamos Scientific Laboratory plasma analyzer, which provides a measure- 
ment of the velocity distribution of solar wind particles with a mini- 
mum time resolution of about 90 seconds, and the Goddard Space Flight 
Center (GGFC) magnetometer, which measures the magnetic field magnitude 
and direction vrlth a time resolution of 6o milliseconds. In addition 
to g, and •;(^,n) other upstream parameters used are the 

upstream niagnetic field magnitude |1^|, the raagnetosonic Mach number M , 
solar wind velocity Alfven speed the sound speed C^, and the 

jjarticle density H. The r.m.s. field strengths for the electromagnetic 
and electros catic components in 5>i2 -seconds average shock spectra like 
those in Figures 1 and 2 are used as the measure of the intensity of 
shock turbulence. We define E and E o r.m.s. electric 

jTins ^ X rnifc3 , d 

field strengths obtained from the average shock electric field spectrum 
by integrating from 20 Hz to 200 Hz, and from 200 Hz to 4 kHz, respec- 
tively. and r.m.s. electric field strengths 

of the electromagnetic and electrostatic components, respectively, of 
the shock electric field spectrum. is the r.m.s. magnetic field 

strength of the electromagnetic component, obtained by integrating the 
average shock magnetic field spectrum from 20 Hz to 4 kHz, The approach 
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of this study is to seek correlation between the values of E 

rms,l^ 

^rms,2^ ^rms^ P. C^, N, 

'!'(^^n)^ and |^j . 

Simultaneous measurement of the spectrum of bow shock turbu- 
lence and upstream solar wind parameters is not possible with a single 
satellite, so the solar mnd parameters derived from the plasma analy- 
zer measurements are averages over time periods during which IMP -6 
ms in the upstream region near the bovj- shock. The solar wind para- 
meters are averaged over one to two hour periods, before or after the 
shock crossing at which the corresponding plasma wave spectrum is 
measured, to obtain values which characterize the solar wind properties 
near the uime of the shock crossing. The one to two hour averaging 
periods also makes the solar wind parameters relatively independent 
of short period fluctuations. If multiple shock crossings take place 
in a time interval less than the averaging times of solar wind para- 
meters, the corresponding series of r.m.s. field strengths for shock 
turbulence is averaged over the number of crossings. The values of 
shock normal angles ^;(^,n) are computed from the model of Fairfield 
[1971] using U-seconds averages of the GSFC magnetometer measurements 
near the leading edge of each shock crossing. 

The distributions of values for the upstream parameters used 
in this study are shown in the histograms of Figures 3 and 4 . The 
general shapes of the distributions compare favorably mth the corre- 
sponding distributions which are obtained for the solar wind when 
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many more samples are taken [Formisano and Moreno, 1974] . The solar 
wand parameters of the present study thus approximate the typical 
ranges of solar wind conditions. The minimums, maximums, means, and 
standard deviations for the distributions of Figures 3 and 4 are 
listed in Table 1, 

Energy densities for the electromagnetic and electrostatic 
components of the bow shock spectrum may be computed from 

\ms,i^ ~ where Gaussian 

units are used for the electric and magnetic fields. The relationship 

^^®rms,i^ solar wind energy density 

e(SW) =|fc^V^^ + |Kg(I^ , 


where m^ is the proton mass and is Boltzmann’s constant, is indi- 
cated by the distributions of the ratios e(E .)/e(SW) and 

shown in Figure 5. It is evident that the plasma wave 
energy densities in the bow shock are al-ways very small fractions of 
the solar vand energy density. The absolute energy density distribu- 
tions for the bovr shock electromagnetic and electrostatic components, 
and for the solar wind are shown in Figure 6. Table 2 list the char- 
acteristic parameters for the distributions of Figures 5 and 6. 

C. Linear Correlation Coefficients 
A standard statistical correlation technique (least squares 
fitting) has been used to calculate linear correlation coefficients 
that indicate the degree of association between the r.m.s. field 
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strength of shock turbulence and solar nind parameters. The fitting 
equation used is log y - bx t a, where y is the r.m.s. field strength, 
X Is the plasm parameter, and b and a are constants. The choice of 
fitting equation is not based on any theoretical relationship between 
the parameters, but is used only to proyide a quantitative measure 
(the correlation coefficient) with which to Judge the statistical 
dependence (or independence) of the parameters involved. Since the 


measurements of electric and magnetic field amplitudes are obtained 
from antenna voltages which are digitised with equal quantising steps 
on a logrithmio scale, the relative uncertainty iv/v in the measured 
signal voltages is constant. Therefore, the least squares fitting is 

done with constant relative unee-r-tain+vr n-v, 

uncertainty in the r.m.s. fiel(3 amplitudes 

[Bevington, 1969, p. 180]. 


Table 3 lists the linear correlation coefficients K obtained 
for the two-parameter fits. For S points used in the least squares 
fit, is the critical value of the correlation coefficient at the If, 
level of Significance for a two-parameter fit [Neville and Kennedy, 
19 S 4 , p. 314 ]. If the absolute value of the computed correlation 
coefficient |H| exceeds E^, the probability is Xf. that the observed 
correlation between the two parameters is due to chanee alone. There- 
fore, a strong correlation is indicated hy a large value of E 

to R . 
c 


compared. 
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III, ELECTROSTATIC TURBULENCE 

A study of the coefficients in Xahle 3 reveals that the electro- 
static field strength is strongly correlated mth several of 

the sclar wind parameters. In particular, 2 and T /T have the 

largest value of R, indicating that a strong positive correlation 
exists between the intensity of electrostatic turbulence in the bow 
shock and the electron to proton temperature ratio in the upstream 
solar wind. Figure 7 shows the plot of E^^ ^ against Ig/T^. The 

diagonal dashed line in Figure 7 is the line of regression for E 

rms,2 

on T^/t^, and the slope of this line indicates the apparent dependence 

^ ^rins,2 large error bars are at + o(y), the standard 

deviation, above and below the line of regression and indicate the 

degree of dispersion of the data points about the mean. However, in 

this case the large dispersion does not mean that large measurement 

errors are present in the values of dispersion of the 

values of is probably indicative of the dependence of E 

’ rms,2 

on other parameters besides T^/t^ which are not included in the two- 

parameter fit. As is shown in Table 3, displays significant 

correlation (R > R^) vrLth other shock parameters, so that the total 

dependence of E^jj,g g is very likely a complex function of T /t , T , C , 

' e' p^ p-’ 

etc. Also, the one to two hours averaging period for the solar wind 
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parameters removes short period fluctuations from these upstream values 

which may he related to some of the scatter in E In addition, the 

rms,2 ^ 

electric field measurements are made with a single antenna so that only 

one component of the electric field is detected at any instant. However^ 

as shewn in paper I, the electric field turbulence in the bow shock 

has a broad angular distribution, and a measurement averaged over 5.12 

seconds (~ half the spin period) can be assumed to be close to the total 

field magnitude. Considering then the approximations involved in a two- 

parameter fit, it can be expected that a true, but weak, dependence of 

^rms,2 °° upstream parameter may not emerge from the correlations in 

this study. However, a strong, and therefore more important, dependence 

oi" 2 ™ upstream parameter should be evident in the correlations, 

even though the functional relationship between E _ and the upstream 

x*nis y ^ 

parameter is not known. The slope of the regression line is the 

important measure of association between E ^ and T /T , and since 

rms,2 e' p' 

the regression line is a mean fit, the standard deviation of the mean, 
a(y), is a better measure of dispersion because this tends to compen- 
sate for correlations not included in a two-parameter fit . The .'^Ttia 1 1 p-r 
error bars in Figure 7 are at the 1^ level of significance, +2.576 o(y). 
At their location on the regression line, the small error bars also 
indicate the 1% l imi ts to the dispersion of the slope when the regression 
line is rotated about its centroid. It is clear that the slope of the 
regression line changes neither in sign nor greatly in magnitude at 
these two limits. 
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The second strongest correlation of E is with the proton 

temperature T^. The plot of E^ ^^ 2 against T^ is shown in Figure 8 , 
where the line of regression for the two-parameter fit is the diagonal 
dashed line and the error bars have the same meaning as in Figure 7. 

The slope of the regression line indicates that a negative correlation 
exists between electrostatic turbulence in the bow shock and the 
upstream proton temperature. Since no significant correlation is found 
between electron temperature T^, it appears that in the 

relationship between 2 the ratio T^/t^ of Figure 7, varia- 
tions in the upstream proton teraparature T is the significant factor. 

tr 

Therefore we conclude that the instability mechanism that produces 
electrostatic turbulence in the bow shock is primarily associated with 
the value of T^/t^, and that the proton temperature T^ serves to modu- 
late the efficiency of that mechanism. Table 3 also shows that E 

rms ,2 

has a significant negative correlation with the sound speed C^, which 

is already implicit in the correlations with T /T and T if we note 

I , ® P P 

that Cg O' [T^ + Tp] • = [T^(T^/Tp + 1)]^' It may also be observed 

that a negative correlation with only marginal 

correlations with other parameters such as and perhaps P, and no 
correlation with i!j(^,n). 

In the column under E^^^ the r.m.s. electric field strength 
between 20 Hz and 200 Hz^ it can be seen that -j^ has its strongest 

correlation with T^/T^; the corresponding points are plotted in Figure 
9 . A comparison of the correlation coefficients in Table 3 reveals 
that E - and E have like correlations with T /T , T , and C , 
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i.e., fath the same sign and ordering. However, . also correlates 

^ JQo j X 

with ^(?,n), unlike o, but instead like B . The correlations 

* X XUS • 

^rms,l thus indicate that the 

dependence of on the upstream parameters includes charac- 

teristics of both purely electrostatic turbulence and purely electro- 
magnetic turbulence. We may therefore conclude that in the lotr fre- 
quency portion of the shock spectrum the electric field is derived 
from the fields of coupled electrostatic and electromagnetic waves . 

A similar conclusion was reached in paper I where it was shown that 
below about 200 Hz the ratio of electric field energy density to 
magnetic field energy density Sg/sg is greater than expected for 
whistler waves, and that therefore the electric field energy density 
eg must also include electrostatic waves. It can thus be seen that the 
entire shock electric field spectrum is primarily composed of electro- 
static waves, which couple to whistlers at frequencies below the elec- 
tron gyrofrequency. An electrostatic wave which can couple to 
whistlers in a high p plasma such as the solar wind is the ion sound 
wave [Formisano and Kennel, I969] . From the discussion of Formisano 
and Kennel [1969] it may be inferred that whistlers and ion sound 
waves can couple over a broad range of frequencies below the electron 
gyrofrequency whenever strong temperature gradients occur; this result 
IS consistent with the observations of electrostatic and electromag- 
netic vraves beloT-j’ about 200 Hz in the bow shock spectrum. 

The observed correlations of and with T^/t^ and 

Tp imply that the kinetic properties of the solar wind plasma, i.e., 
the wave-particle interactions that modify the electron and proton 
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velocity distributions, are the primary factors determining the in- 
tensity of electrostatic turbulence in the bow shock. Usually T > T 

e p 

in the upstream region and leading edge of the shock. The observation 

of temperature jumps in T and T across the shock such that T is 

P p 

2-4 times greater than in the near downstream region means that pro- 
tons provide the major dissipation mechanism for the bovr shock [Montgomery 
£4 1970] . The electric field spectrum from which EL , and 

4l*lDS ^ J- 

^rms,2 computed can thus be associated with jaroton heating. 

Such fluid parameters as M^, p, and KB,n), with which E 

" T 7 j rms,l 

and weak or null correlations, are apparently not impor- 

tant determinants of the intensity variations in the shock electric 
field spectrum. However, it should be noted that the m.h.d. structure 
of the bow shock is often classified in terms of M^, p, and ,i,(B,n) 

[Pormxsano and Hedgecock, 1973a; Greenstadt, 1974], and within these 
classifications, almost all of the shocks used in the correlations of Table 3 
fall into the category in which the bow shock is most often found, that 
of "turbulent shocks" with > 3, P > 0.1, and i!r(^,n) ~ 30“ - 70°. 

Certain critical values of the upstream fluid parameters, M* ~ 2.5 - 
3.0, p ~ 0.1, and il((B,n) = 88.7“, at which significant changes in the 
m.h.d. structure of collisionless shocks occur [Woods, 1969a, b, 

I97O; E^ul, 1972; Tidman and Krall, 1971j chap. 3j Biskamp, 1973] 
not well represented in the correlations of Table 3. 
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IV. ELECTROMAGNETIC TURBULENCE 

For the upstream parameters of Table 3> electromagnetic tur- 
bulence, represented by the magnetic field strength B , is found 
to display a strong correlation v/ith the particle density R. Figure 10 
shows the plot of and N; the line of regression and error bars 

have the same meaning as in previous figures. The slope of the re- 
gression line indicates that more intense magnetic field fluctuations 
occur for larger values of the particle density, a variation consistent with 
whistler waves since the index of refraction (~ cB/e) for the whistler 
mode increases with particle density. The positive correlation 
between B^^^ and N supports the previous conclusion that whistler 
waves comprise the shock magnetic field spectrum. The slope of the 
regression line in Figure 10 is not drastically changed if the relatively 
few points at high density are omitted from the fit . 

The shock normal angle \];(B,n) is the only other parameter of 
Table 3 mth which B^^^ shoves a substantially significant correlation. 

The plot of and i!fC5,n) is shown in Figure 11. 'f»(^,n) takes values 
in the range 0° to 90 % with ^(2,n) = 0° being defined as a parallel 
shock, 0“ < < 88.7“ is an oblique shock, and 88.7“ < '!iCS,n) ^ 

90“ is a perpendicular shock. These definitions for iii(S,n) are 
theoretical [Tidman and Krall, I 97 I, p. 24]; the experimental values 
for shock normal angles have error limits of + (5° - 10“) so that 
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most shocks can he considered as oblique, ~ o"to 85% or per- 
pendicular, ~ 85“ to 90°. Other classifications are also 

used [Greenstadt, 1974] . Whistlers propagate in a cone about the 
direction of % therefore as t lies closer to the plane of the shock, 
i.e., as tj((3,n) 90% then on.the average it can be expected that 

the energyr density of whistler .^ves in the plane of the shock will 
increase. Therefore, the positive correlation of and '^(S,n) 
also supports the identification of as whistler turbulence. 

The frequency range of B^^ (20 Hz - 4 kHz) corresponds to . , 
high frequency end of a much broader spectrum of magnetic field 
fluctuations of the bow shock region which generally has its largest 
spectral densities at frequencies near and below the proton gyro- 
frequency (0.001 - 1.0 Hz) [Holzer et al., I966, I972; Olson et al., 
1969; Fairfield and Hess, I97O] . Since most of the magnetic field 
energy density is at the much lower frequencies associated with the 
m.h.d. regime, it is not surprising that B^^^^ has insignificant 

correlations with m.h.d. fluid parameters such as M and p. 

A. 

A positive correlation between B^^^ and the upstream magnetic 
field magnitude |bJ is also indicated in Table 3. The observed corre- 
lation is not related to the shift of the Mistier cutoff, i.e., the 
electron gyrofrequency , toward a higher frequency in the spectrum of 

magnetic field fluctuations since for a typical spectrum that varies 
-4 

as f the largest expected upward shift in the whistler cutoff increases 

®rms about 1^, whereas the values observed for B range over 

3rins 
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more than an order of magnitude. The correlation of B with IbI 

rms ' I 

thus indicates an actual statistical increase in the intensity of 
magnetic field fluctuations with an increase in jBj . The correlation 
appears to he weak, however. 
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V. BOW SHOCK STRUCODRE 

The leardng edge of a coUisionless plasma shock -wave has been 
considered in theoretical studies to he a region vSiere particle re- 
flection and heating can occur Tjhich greatly influence the dissipation 
processes further into the shock [Woods, 1969b, 197I; Biskamp, 1970, 

1973; Biskamp and Welter, 1972; Tidman and, Krall, I971, p. I30] . 
Experimental measurements of electron and ion velocity distributions 
near the earth's bow shock have shown that electron preheating and 
ion deceleration does occur near the foot of the magnetic field 
gradient [Montgomery et al,, 1970; Weugebauer, I97O; Weugebauer ^ al., 
19713 . The leading edge and transition region of the bow shock can often 
be resolved in the structure of the plasma wave spectra that are ob- 
tained at a given shock crossing. It may therefore be ejcpected that 
electron preheating should result in a broad spectrum of turbulent elec- 
trostatic waves at the leading edge of the bow shock. In fact, electro- 
static mves are typically observed on the upstream side of the transi- 
tion region in almost all bow shock crossings that we have studied. 

In the main transition region of the shock the spectrum of plasma 
waves is modified significantly and can be associated with the ion 
thermalization that provides the major dissipation for the shock 
[Montgomery et al. ., 1970] . 
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The data in Fif jre 12 were obtained for a perpendicular shock 
with the transition region clearly defined. The upper panel shows 
the electric field in the 3. 11 kHz channel of the University of Iowa 
spectrum analyzer which was measured ^rith the rapid-sample mode 
described in paper I. The magnetic field measured by the GSFC magneto- 
meter is shown in the lower panel. The angles and 0^ are the solar 
ecliptic longitude and latitude of 5. Before about l440:38 UT the 
spacecraft is in the quiet upstream solar mnd. At l440:38 UT the 
leading edge of the bow shock is encountered as shown by the sudden 
increase in electric field strength above the solar wind noise level 
and the corresponding onset of fluctuations in the magnetic field. 

It can be seen that the electric field increases substantially in 
the time interval between the vertical dashed lines @ and , 
beiore the main gradient of the magnetic field which occurs in the 
time interval between vertical dashed lines 0 and 0 . Since 
the spectrum of electron plasma oscillations (peaked at about f ~ 

30 kHz) IS often observed to broaden toward 3-11 kHz upstream of the 
bow shock, the electrostatic noise between points (0) and (0^ in 
Figure 12 may ini5.eate a heating of the solar wind electron distribution 
just before the main gradient of the magnetic field oecurs. 


Such preheating of electrons would substantially increase the value 
of Tg/Tp in the leading edge of the shock with a consequent lowering 
of the threshold for the destabilization of drift currents such as 
occur with ion sound waves [Tidman and Kfall, I971, p. 128], The 
electrostatic turbulence between (b) and (0 correlates well with the 
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magnetic field gradient, similar to the resu 3 .ts of Fredricks et 
^ [ 1968 , 1970 b] . It can be ejqpected that most of the ion thermaliza- 
tion must occur in the region of this large gradient [Montgomery ^ al., 
I 97 OJ Greenstadt ^ a-1. , 1975] • It is clear that in the shock cross- 
Figure 12, the structure of electrostatic turbulence in the 
leading edge of the shock is not greatly different from that in the 
main transition indicating that the regions of electron and ion 
heating overlap. In the immediate downstream region, the electric 
field noise intensity at 3-H kHz shows a series of nearly periodic 
fluctuations that display exponentially decreasing magnitude. The 
periodicity of the larger downstream fluctuations is not related to 
the spacecraft spin period, so it may be assumed that the fluctuations 
are real time variations of electrostatic turbulence. The exponential 
envelope of the downstream fluctuations implies that the spectrum of 
electrostatic waves is being damped, perhaps by the type of fluctua- 
tions in the ion distribution observed by Montgomery et al. [ 1970 ] 
in the downstream region. 

A second example of a perpendicular shock in which the structure 
of the transition region is clearly defined is shown in Figure 13 . 

The upper panel shows the electric field in four channels (36 Hz, 

311 Hz, 3*11 kHz, 31-1 kHz) of the plasma wave spectrum analyzer, 
with average measurements plotted as vertical bars and peak measurements 
indicated by dots. The magnitude and direction angles of the 
magnetic field measured simultaneously by the GSFC magnetometer is 
shown in the lower panel. A vertical dashed line at about 2219:00 UT 
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marks the time at which the shock transition region begins as determined 
by the electric field, measurements. It can be seen that the electric 
field noise begins almost one minute before the main gradient of 
the magnetic field. The shock transition region is fairly wide in 
time and the 5.12-seconds average measurements can resolve the changes 
in the electric field spectrum as the transition region is traversed. 

The beginning of the transition region is detected simultaneously in 
all channels; however, the rise times to maximum electric field in- 
tensity in the transition region is progressively shorter for increasing 
frequency. The relaxation from maximum intensity in the two high 
frequency channels begins before the maximum intensity is reached 
in the two low frequency channels so that by the time the maximum field 
strength occurs at 3U Hz, the electric field intensity at 3I.I kHz is 

nearly back to the noise level. The magnetic field measurements 
indicate that the electrostatic turbulence at lower frequencies cor- 
relates with the main transition region vh±la the electrostatic tur- 
bulence at high frequencies correlates with the leading edge of the 
bow shock. 

Between 2219:20 and 2220:00 UT, when the spacecraft was in the 
leading ..dge of the shock transition, a rapid-sample measurement of 
the 31.1 kHz electric field channel occurred. The rapid-sample data 
is shown in the polar plot of Figure l4 in which the electric field 
measurements are plotted versus the spin angle of the sjacecraft 
antenna in the solar ecliptic plane. During the time of the rapid- 
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sample measurements the average magnetic field direction in the lead- 
ing edge of the how shock had a solar ecliptic latitude 0 ~ 30° and 
a solar ecliptic longitude ~ 90°. The projection of the magnetic 
field vector into the solar ecliptic plane at ~ 90° is indicated 
in the polar plot of Figure l4. By computing the second moments of 
the rapid-sample measurements about the solar ecliptic x and y axes 
(equivalent to evaluating the moments of Inertia for unit "mass" 
particles with moment arms equal to the fluctuation amplitudes) and 
diagonalizing the resulting matrix, the principal axes may he obtained 
for the two-dimensional distribution of rapid-sample measurements. 

The direction of the major principal axis, indicated by on the 
polar plot, represents the average polarization direction of the 
electric field fluctuations. It is evident that the average electric 
field polarization is nearly aligned with the magnetic field vector 
directxon (0^ - ~ 10°); the frequency and average polarization 

therefore identify the noise as electron plasma oscillations [Fredricks 
et 1968, 1970b; Rodriguez and Gurnett, 1975]. The spectrum of 
electron plasma oscillations at about 2219:00 UT (20 seconds before 
the rapid-sample data of Figure 14 was obtained) is sharply peaked at 
the electron plasma frequency near I6.5 kHz. The broadening of the 
spectrum of electron plasma oscillations toward the higher frequency 
31.1 kHz indicates that heating of the electron distribution began 
just before the rapid-sample meastirement was taken, and occurs throu^- 
out the rapid-sample time interval. The five consecutive 


average 
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electric field spectra measured between 2220:00 and 2220:26 UT imme- 
diately after the rapid-sample measurement of the 31.1 KHz channel and 
corresponding with the main gradient of the magnetic field (see Figure 
13) are shown in the lower panel of Figure l4. Each spectrum is 
labeled by a time-ordered snapshot number. The averaging time of 
snapshot 1 includes the foot of the magnetic field gradient. Snapshot 
5 is the most intense electric field •'’pectrum of the sequence and is 
defined as the shocfc electric field spectrum, similar to those sho^ 
in Figure 1. The development of the broad peaK centered near the 
3 oil JsHz channel in snapshot 1 correlates with the appearance, in 
the leading edge of the shock, of 0.5 - 4.0 Hz whistler waves of the 
type reported by Fairfield [1974] . Further upstream, vdiere the 
spectrum of electron plasma oscillations is sharply peaJted at the elec- 
tron plasma frequency near I 6.5 kHz, the 0.5 - 4.0 Hz waves have damped 
out. These upstream whistler waves probably involve resonant inter- 
actions with electrons [Fairfield, 1974] . Therefore the broadening of 
the spectrum of electron plasma oscillations toward lower frequencies 
probably indicates scattering interactions between upstream electrons 
and the 0.5 - 4.0 Hz whistler waves. The overall broadening of the 
spectrum of electron plasma oscillations toward 3 I.I kHz and 3 .II kHz 
thus clearly indicates the occurrence of electron heating in the leading 
edge of the shock and we interpret the peak at 3 .II kHz in snapshot 1 
to mean that maximum heating of the electron distribution has occurred. 
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Following the sequence of snapshot ntanbers 1-5, it can he seen 
that as the min gradient of the mgnebij^iilLeld of the shock is 
traversed the electric field spectrum develops a peak at 311 Hz simul- 
taneously with the disappearance of the peak at 3-11 kHz. Since proton 
heating in the how shock is observed to occur after electron heating 
[Montgomery et al., 1970; Formisano and Hedgecock, 1973h], i.e., 
deeper within the shock transition region, we therefore conclude that 
the broad peak centered at 311 Hz in the shock spectrum corresponds to 
the occurrence of maximum proton heating. This broad peak of electro- 
static turbulence in the shock spectrum is usually centered betvreen 
200 - 800 Hz, corresponds to ^ in the correlations of Table 3, 
and is characteristic of almost all shock crossings, as shown by the 
sample Ox shock spectra in Figure 1. Also characteristic of most bow 
shock crossings is the sequential development in electric field spectra 
through the shock transition shown by the snapshots 1-5 of Figure Ik, 
which suggests that the electrostatic turbulence associated with the 
3.11 kHz peak couples into lower frequencies during the proton heating 
process. It is interesting to note that this variation in electro- 
static turbulence with frequency indicates that a strong coupling 
mechanism exists between the electron and proton heating processes 
in the bow shock since the direction of energy cascade in the spectrum 
of electrostatic turbulence is toward lower frequencies rather than 
toward higher frequencies as mi^t be expected [Hoth, 1971] in a 
turbulent spectrum analogous to conventional hydrodynamic turbulence. 
Ihst the downstream edge of the transition 


region the uppei* panel of 
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X3 sliows *fcli8#"b 'c2i0 GXsc*fa3T os'fcs.'fcic "fcixrXiiXGiics B.*fc Xow fjT0{^u.0Xicx€s 
decareases Taut stiXX remains at an intensity that is Xarge compared, 
with the noise at higher frequencies. The resuXting eXectric field, 
spectrum in the near downstream region resemhXes the shock spectrum 
at a reduced intensity. 


It can aXso be seen in Figure X3 that in the upstream region 
before about 2213 UT the GSFC magnetometer data shows the presence of 
long period (3-60 seconds) waves of the type which are thought to be 
generated by superthermaX proton streams reflected from the bow- 
shock [Greenstadt efc aX. j X968, X97O; Fairfield, 1969; Barnes, 1970; 
Scarf X97O; Fredricks, 1975]. Scarf [1970] have 

observed direct correlations between electrostatic oscillations at 


3 Mfe, associated with the ion plasisa frequency in the solar vjind, 
and the long period magnetic field oscillations . The top panel of 
Figure . 13 shows that electrostatic oscillations at 31.I kHz, associated 
with the electron plasma frequency, are also in direct correlation with 
the long period -tTaves in the interplanetary magnetic field. The 


electrostatic oscillations at 31«1 kHz prohahly result from the electron 
stream that must accompany the reflected proton stream to provide 
charge neutralization. 
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VI. SUMMAEY OF EESDIIS Alffi CONCLUSIONS 


THe plasm ^ve turbulence of the earth's bow shock displays 

distinct correlations with several upstream solar wind, parameters. 

It is found that the intensity variations of electrostatic turbulence 

in the bow shock corredates strongly with the ratio of electron to 

proton teii^eratures and with the proton temperature T , as 

measured in the upstream solar wind. For large values of T /t , 

large values of the electrostatic r.m.s. field strength E occur 

rms,2 ^ 

where is obtained by integrating the shock electric field 

spectrum from 200 Hz to 4 kHz. The negative correlation found between 
^rms,2 ^p^ implies that changes in the upstream proton tempera- 
ture modulate the efficiency of the electrostatic turbulence mechanism 
which is associated with the value of if we make the plausible 

assumption that more intense levels of electrostatic turbulence occur 
Tdien an unstable plasma is further removed from the tlxreshold of 
stability, then we can relate the observed correlations of 2 with 
^e'^^p ^p specific electrostatic wave modes whose stability 
criteria depend on T^/t^ and T^. Two-stream instability criteria 
are often expressed in terms of the two parameters and v^/v , 

where v^ is a relative drift velocity and v^ is the electron thermal 
speed [Stringer, 1964 ; Tidman and Krall, I97I, chap. 7]. A discussion 
of two-stream instability modes that are candidates for bow shock 
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turbulent heating ineehanisms has been given by Greenstadt and 

Fredricks [1974], who derive stability criteria in terms of the 

minimura relative drift speed v^ between electrons and ions which must 

be exceeded over the scale length of the shock magnetic field gradient 

AB/b in order for instability to occur. Various scale lengths are 
o 

chosen, all of which are proportional to and/or (- 8nMni^/}B[ ). 

The temperature ratio T^/T^ and the relative drift velocity v^ (or 

respectively indicate the relative widths of the electron and 
ion distribution functions and the relative displacement between 
the nw.yimiinif! of the two distributions . In terms of and v^ the 

two-stream stability criteria indicate how well resolved from each 
other in velocity space the two streams have to be before unstable 
waves are generated. In general, anything that increases the resolu- 
tion between the streams enha.nces the potential for instability and 
would be expected to lead to more intense levels of electrostatic 
turbulence. In particular, the positive correlation between 2 

and T /t suggests that electrostatic turbulence in the bow shock is 
e p 

associated with a two-stream instability, the stability threshold of 
which is probably exceeded if T^ increases substantially in the leading 
edge of the bow shock. Electron preheating in the leading edge is 
confirmed by the broadening of the spectrum of electron plasma oscil- 
lations over the range of frequencies 3-30 kHz. For typical solar wind 
conditions only the electron distribution undergoes an appreciable 
drift in the bow shock gradient [Wu and Fredricks, 1972] so that the 
instability arises from electrons drifting through an ion background. 
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Tlis 37fiS'uX*b of iiliis siiudv ijliai* 7? _ * 4. 

^ rms,2 ^ strong negative correlation 

with the proton temperature is consistent mth an electron-proton 
streaming instability because for a given relative displacement v^/v^ 
a smaller vidth to the proton distribution f^Cv) contributes to babtL 
resolution of the two streams and an enhanced potential for instability. 
Hius, -variations in could serve to modulate the intensity of 
the resulting electrostatic turbulence. 

As indicated by Greenstadt and Itedricks [1974], streaming 
betw.cn electron and proton distributions preferentially heats electrons, 
which corresponds to the observation of aaidmnm electron heating 
occurring in the leading edge of the bow shock. To heat ions, the 
major dissipation mechanism of the bow shock, to temperatures such 
that Tp ~ (2-4)Tg requires a streaming between ions [Auer ^ al., I971; 
Ihpadopoulos, 1971; Biskamp and Welter, 1972], or other nonlinear 
instability modes. In the correlations of Table 3, corresponds 

to the electrostatic turbulence observed in the main shock gradient 
and was identified with the maximum heating of protons. The strong 
correlation of with and T^ as measured upstream of the 

bow shocic thus implies that the heating of protons in the main transi- 
tion is strongly coupled to the heating of electrons in the leading 
edge of the shock. The coupling of the electron and proton heating 
mechanisms is most clear]y shown by the characteristic sequential 
development of the electric field spectrum through the shock tiansi- 
tion, an example of which is shown in Figure l4, in which the intensity 
of electrostatic turbulence at lower frequencies 300 Hz) increases 
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simultaneously vath a decrease in electrostatic turbulence at the 
higher frequencies associated with electron preheating. The low 
frequency (200-800 Hz) peak in the spectrum of electrostatic turbulence, 
associated with the maximum in proton heating, is characteristic of 
almost all shock electric field spectra. 

low frequency (20 Hz - 200 Hz) electric field turbulence, 

represented by is found to correlate strongly with and 

Tp, similar to the correlations of with T^/T^ and T^, and 

therefore indicating that a substantial portion of the electric field 

energy- density at low frequencies is derived from electrostatic waves. 

The entire electric field spectrum measured in the bow shock is thus 

primarily composed of electrostatic turbulence, and is associated 

with the occurrence of maximum proton heating. Since E also 

rms, 1 

correlates with the shock normal angle ^(B,n), similar to the correla- 
tion of mth i!f(B,n), this is taken as evidence that E is also 

dei-ivea from the eleotrio field of an eleotrcmghetia mode, probably 

vhiatlera, which couplea to electrostatic wares at frequcacies below 
the electron gyrofrequency . 

Electromagnetic turbulence in the bow shock, represented by 
the r.m.s. field strength of magnetic field fluctuations B in the 
range 20 Hz - 4 kHz, is found to show positive correlations with the 
particle density W and the shock normal angle KB,n). These correla- 
tions are consistent with whistler turbulence with high density, per- 
pendicular shocks having the lai^gest values of Upstream fluid 

parameters, such as and P, used in the m.,h.d. description of the 
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ho-^ shock do not correlate vith ireanmbly because the range 

of frequencies for covers only the high frequency, low spectral 

densila' portion of the total Mgnetic field spectrum near the bow 
shoek. 
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Table 1 

Solar Wind Ibraaeter Characteristics 


Min Max Mean Standard 

Deviation 



1.5 

26.6 

6.8 

4.7 

3 

0.03 

4.5 

1.02 

0 . 8 ‘i 

|s| ( y) 

2.0 

16-9 

7.5 

3-2 

^s 

3.9 

12.5 

7.8 

1.9 

Sg (ion/s) 

33.2 

105 . 

58. 

14 . 

W (cm 

0.5 

19.0 

5.2 

3.4 


15.4 

278. 

88. 

60. 

t(B,n) (deg) 

11.7 

90. 

60. 

19. 

T /t 

e P 

0.65 

9.0 

2.6 

1.8 

® ) 

315 . 

660. 

431. 

63. 

Tg (x 10^ °K) 

0.6 

6.0 

1.6 

0.80 

(x 10^ °K) 

Jr 

0.15 

3.0 

0.98 

0.72 
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Table 2 

Energy Density Characteristics 


Min Max Mean* Standard 

Deviations 


(x 10"®) 

0.0028 

6.4 

0.17 

0.021 

1.3 

'(®rn.a,2Ve(SW) 

(X 10"®) 

0.0001 

14.2 

0.054 

0.0047 

0.63 

(x 10"^) 

0.026 

5.9 

0-38 

0,10 

1.4 

e(E .) 

(x 10 ^ergs cm"^) 

3.7 

3020. 

105. 

16. 

672. 

®^®rms_,2^ 

(x 10 ergs cm"^) 

0.32 

6120. 

34. 

2.9 

4oo. 

e(B ) 
rms'^ 

(x 10 ergs cm"^) 

1.1 

2600. 

24. 

5.0 

115. 

e(SW) 

(x 10 ergs cm"^) 

6.9 

44o. 

59. 

22. 

161. 


•X- 

Logarithmic means and standard deviations 


0 


-0.300 

-0.266 


-0.294 

-0.286 


"a 


0.286 



0.342 

l^i 


0.306 

M 

s 

0.186 


Te/Tp 

0.4l6 


V 

sw 

-0.165 


C 

s 

-0.328 


N 

0.032 


Te 

-0.200 


Tp 

-0.376 


S 

81 

96 

RcW 

0.281 

0.259 


0.236 

0.599 

-0.361 

-0.463 

0.242 


-0 .241 
-0.566 
75 


0.136 
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Figure 1 


Figure 2 


Figure 3 


FIGURE CAPTIONS 



The intensity variability of bow shock electric 'field 
turbulence is indicated by the representative sample of 
shock spectra obtained in 36 crossings of the bow shock. 
The average spectra are 5.12 seconds averages and the 
peak spectra are 0.1 second aveiages. E is the r.m.s 
electric field strength obtained by integrating a given 
spectrum from 20 Hz to 200 kHz. 


The iirbensity variability of bow shock magnetic field 
turbulence is represented by average and peak spectra 
obtained at the same 36 shock crossings of Figure 1. 

®rms r.m.s. magnetic field strength. 

The distributions of values for solar wind parameters 
used in the present study. The mean and standard devia- 
tions for each distribution are indicated by the dot and 
error bars, respectively. 

The distributions of values for solar wind parameters used 
in the present study, similar to Figure 3. Means and 
standard deviations indicated as in Figure 3. 


Figmre 4 


Figure 5 


Figure 6 


Figure 7 


lj-5 

The distributions of the ratios of plasma wave energy 
densities to the solar wind energy density. The logari- 
thmic mean values and standard deviations of the ratios 
are indicated on the plots . 

The distributions of the absolute values of plasma wave 
energy densities and solar wind energy density. 

The plot of the r.m.s. field strength of the elec- 

trostatic component of the shock spectrum (200 Hz - 4 kHz), 
against T^/t^ which shows a strong positive correlation. 

The dashed diagonal line is the line of regression for a 
least squares fit to the equation log y = bx + a, where 

^ ®rms,2 ^ “ ^e^^p' slope of the regression 

line, which is the important meausre of association, has a 
dispersion which is indicated by rotating the regression 
line about its centroid to the limits indicated by the 
small error bars at + 2.576 o(y). The large error bars 
at + a(y) probably arise from short period fluctuations 
fi*oni correlations mth other upstream para- 
meters not considered in this two-parameter fit. 
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Figure 8 The plot of 2 the upstream proton tempera- 

ture Tp shors^ing a strong negative correlation for a 

least squares fit to log y = bx + a where v = E 

rms,2 

and X = T^. The error bars have the same meaning as in 
Figure J. 

Figure 9 The plot of against T^/t^ which indicates a posi- 

tive correlation. E^^ ^ is the r.m.s. fi.eld strength 
of the electric field between 20 Ha and 200 Hz in the 
shock spectrum, and is observed to correlate with T /t 

e' p 

similarly to the correlation shown in Figure 7 between 

^rms,2 ^e/^p' error bars have the same meaning 

as in the two previous figures. 

Figure 10 Ihe plot of the r.m.s. magnetic field strength of 

the electromagnetic component of the shock spectrum 
(20 Hz - 4 kHz), against the solar wind particle density 
N. A positive correlation is indicated. The error bars 
have the same meaning as in previous figures. 

Figure 11 The plot of against the shock normal angle ^(B,n), 

showing a positive correlation- The error bars have 
the same meaning as in previous figures. 



4 ? 


Figure 12 


Figure I3 


Eapid-sample measurement of the electric field at 3. 11 
kHz showing the correlation •vrith the transition region 
as determined hy the magnetic field. She interval 
between the vertical dashed lines @ and (g) indicates 
the leading edge of the shock where the initial electro- 
static noise occurs. The main transition region, between 
® and © , corresponds to the large gradient in the 
magnetic field and the associated electrostatic turbulence, 
the intervals © - ® and ® - © probably correspond, 
respectively, to the regions of mayimnni electron and 
proton heating. Downstream fluctuations in electrostatic 
noise display an exponential damping which may result from 
further proton heating. She apparent periodicity of the 
major downstream fluctuations is not related to the 
spacecraft spin period so it may be assumed that the 
fluctuations are real time variations- Shock parameters: 

Ma = ^'8. 3 = 0.28, TjT ^ = 6.7, ,^CB,A) = 87“. 

A shock crossing which shows the relation between the 
rise and relaxation times of electric field turbulence 
at high and low frequencies . Electric field noise is 
clearly detected ahead of the main gradient in the 
magnetic field, ihe solid vertical lines indicate the 
time intervals for the rapid-sample measurement and 
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electrxc field spectra of Figtire l4. Upstream electro- 
static oscillations at 3.11 kSz and 31.1 kHz are 
correlated witH the long period -waves in the magnetic 
field. Shock parameters; = 6.9, g = 1.42 , T /t =r 1.0^ 
,jf(^,n) = 90°. ^ 

Figure l4 The polar plot is a rapid-sample measurement plotted 
versus the spacecraft spin angle in the solar ecliptic 
plane^ and represents the electrostatic turhulence at 
31.1 kHz in the leading edge of the shock crossing in 
Figure 13. The orientation of the upstream magnetic 
field vector is indicated hy The major principle 
axis of the distribution of rapid-sample measurements 
xs xndicated by 0^. The electric field spectra in the 
lower panel indicate the sequential developnent of 
electrostatic turbulence with frequency as the main 
transition is traversed. 
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ELECTRON TO PROTON TEMPERATURE RATIO, Te/,tp 
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ELECTRIC FIELD, 200 Hz - 4 kHz 
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Figure 8 
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